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of a Metal Matrix Composite at Elevated Temperatures 
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Mark E. Tuttle, University of Washington 


ABSTRACT 


This research investigates the response of a fiberless thirteen layer hot isostatically 
pressed Ti-15-3 laminate to creep, constant strain rate, and cyclic constant strain rate 
loading at temperatures ranging from 482°C to 649°C. Creep stresses from 48 to 260 
MPa and strain rates of .0001 to .01 m/m/sec were used. Material parameters for 
three unified constitutive models (Bodner-Partom, Miller, and Walker models) 
were determined for Ti-15-3 from the experimental data. 

Each of the three models was subsequently incorporated into a rule of mixtures and 
evaluted for accuracy and ease of use in predicting the thermoviscoplastic response 
of unidirectional metal matrix composite laminates (both 0° and 90°). The laminates 
were comprised of a Ti-15-3 matrix with 29 volume percent SCS^ fibers. The 
predicted values were compared to experimentally determined creep and constant 
strain rate data. It was found that all three models predicted the viscoplastic response 
of the 0° specimens reasonably well, but seriously underestimated the viscoplastic 
response of the 90° specimens. It is believed that this discrepancy is due to compliant 
and/or weak fiber-matrix interphase. 

In general, it was found that of the three models studied the Bodner-Partom model 
was easiest to implement, primarily because this model does not require the use of 
cyclic constant strain rate tests to determine the material parameters involved. 
However, the version of the Bodner-Partom model used in this study does not 
include back stress as an internal state variable, and hence may not be suitable for 
use with materials which exhibit a pronounced Baushinger effect. The back stress is 
accounted for in both the Walker and Miller models; determination of the material 
parameters associated with the Walker model was somewhat easier than in the 
Miller model. 
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CHAPTER 1 INTRODUCTION 


1.1 Background 

In recent years, metal matrix composites (MMCs) have been widely investigated 
for use in satisfying structural requirements in high temperature, high performance 
aerospace applications. Titanium alloys reinforced with silicon carbide fibers are 
currently being considered for use on the National Aerospace Plane (NASP). This 
type of MMC combines a relatively ductile metal matrix material with stiff, brittle 
ceramic fibers. Although the thermoviscoplastic behavior of titanium and other 
metals has been modeled by a number of researchers 12, 7-9, 11, 14, 16-18, 25, 27, 29-31, 
42-44, 47-48], the presence of stiff, brittle fibers significantly alters the behavior of 
MMC's. To date, little has been done to model the thermoviscoplastic behavior of 
MMC's at elevated temperatures. 

The specific MMC studied in this work has a Ti-15-3 matrix reinforced with 
continuous silicon carbide (SCS^) fibers. Ti-15-3, a metastable beta strip titanium 
alloy, is ageable to a high strength and is relatively insensitive to certain corrosive 
environments. The Ti-15-3 matrix undergoes significant inelastic deformation at 
elevated temperatures. In contrast, during this study it will be assumed that the 
SCS6 fibers deform elastically, even at elevated temperatures. 

One of the central issues of this investigation is the selection of a constitutive 
model capable of accurately predicting the thermoviscoplastic behavior of Ti-15-3. 
Three thermoviscoplastic constitutive models— the Bodner-Partom, Miller, and 
Walker models-have been selected for study in this work. The goal of the research 
described herein is to develop a combined experimental-analytic methodology for 
predicting the thermoviscoplastic behavior of unidirectional MMC laminates at 
elevated temperatures. 

Because each of these models has previously been used to describe the behavior of 
one or more thermoviscoplastic materials, it is reasonable to expect that the 
thermoviscoplastic behavior of Ti-15-3 can be adequately predicted by a modified 
form of one or more of these models. Furthermore, the experimental parameters 
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associated with each of these models can be derived from roughly the same database. 
By generating an appropriate database for Ti-15-3, the predictive power of each of the 
models can be assessed for the "neat" (fiberless) material. Each model can then be 
integrated into an appropriate form for a unidirectional laminate using a "rule of 
mixtures" approach. 


1.2 Project Phases 

This study was initiated in February 1989 via a grant (NAG-1 -974) from the NASA 
Langley Research Center (LARC). The project was completed in three phases. 

• Phase 1 : Phase 1 consisted of extensive testing at both UW and LARC using neat 
Ti-15-3 specimens. During this phase an extensive database for the thermo- 
viscoplastic response of Ti-15-3 at elevated temperatures was established, and 
parameters for the three thermoviscoplastic constitutive models were determined. 
Portions of this work have been reported previously by Rogacki and Tuttle [34-35, 
45). 


• Phase 2: During Phase 2 the thermoviscoplastic behavior of unidirectional 
SCSs/Ti-15-3 MMC's was measured at UW and LARC. 


• Phase 3: The "rule of mixtures" formulation (described in Chapter 4) was used to 
predict the thermoviscoplastic response of unidirectional SCS6/Ti-15-3 MMC's 
during Phase 3. This entailed combining the results of Phases 1 and 2 and comparing 
the predictions with the measurements obtained during the MMC tests. 


CHAPTER 2 DESCRIPTION OF LAMINATE MATERIALS 


2.1 Laminate 


The specific MMC studied was manufactured by Textron Specialty Materials of 
Lowell, MA, using a Timet Ti-15-3 matrix (a shortened designation for Ti-15V-3Cr- 
3Al-3Sn) reinforced with continuous SCS^ fibers made by Textron. The neat 
laminates consisted of thirteen layers of Ti-15-3 foil joined by hot isostatic pressing 
(HIP). Although the fiber volume fraction (Vf) was specified as 32%, subsequent 
analysis of micrographs revealed that the Vf was closer to 29%. (See Figure 2.1. Note: 
all figures are located at the end of the respective chapters.) The fibered laminates, 
eight layers thick, had alternating layers of Ti-15-3 foil and continuous SCSg fibers 
joined by hot isostatic pressing. 

For the case of strong, brittle SCS^ fibers in a more ductile Ti-15-3 matrix, the fibers 
will deform elastically while the matrix deforms first elastically then viscoplastically. 
As long as the fibers are not degraded by a chemical reaction with the matrix and the 
composite contains more than a certain minimum volume percent of fibers, the 
composite will fail at a total strain approximately equal to the strain of the fibers at 
their ultimate tensile strength [40]. Therefore, the fibers are treated as being purely 
elastic over the strain ranges encountered in this analysis, and the viscoplastic 
behavior of the MMC is assumed to be directly attributable to the matrix material. 


22 Matrix 

Ti-15-3 is a metastable beta strip titanium alloy. The alloy is weldable, ageable to 
high strength, and relatively insensitive to certain corrosive environments. Ti-15-3 
also offers cost advantages when compared to some other titanium alloys with 
comparable properties [38]. 



Ti-15-3, like other metastable (3-Ti alloys, was developed to satisfy the need for strip 
producible, cold formable, high strength titanium alloys [4J. p-Ti alloys have 
suppressed P transus temperatures and can be quenched from above the material's 
P-transus without martensitically decomposing the P-phase. Metastable P-Ti alloys 
decompose into an a+P mixture upon aging, wherein the a phase has a hexagonal 
close packed crystalline structure and the p is body centered cubic. Ti-15-3 is solute 
rich, and, as such, nucleation kinetics are slower than for solute lean alloys. Also, 
solute rich p-Ti alloys do not form the athermal omega phase. In Ti-15-3, the 
addition of tin as an alloying ingredient helps suppress omega formation [4]. Grain 
boundary a is present in essentially all p-Ti alloys, limiting their toughness [15]. 

In order to optimize ductility, Ti-15-3 is delivered in its metastable beta form. 
After shaping is completed, the alloy is strengthened (and simultaneously made 
more brittle) by causing it to a undergo a phase transformation to its stable alpha 
phase. This aging process is accomplished at temperatures 150-350°F below the P 
transus (1400 ± 15°F) for varying lengths of time, depending on the aging 
temperature. In this temperature regime, the a phase precipitates directly [3,20]. 


Figure 2.2a is a photograph of the "as received" neat Ti-15-3 laminate at 200x 
magnification. The large beta grain structure is readily apparent. The material was 
mounted, sanded, polished, and etched with Kroll's reagent prior to the photograph. 
The dark spots within the photograph are most likely small clusters of the alpha 
phase which were precipitated during the HIP process. (Some of the spots may also 
be voids or contaminants within the material.) Since in many places the beta grain 
boundaries extend across the original lamina interface, the specimens seem to have 
been very well consolidated during the HIP process. However, the fact that the 
interfaces are readily apparent after etching indicates preferential alpha formation 
and/or the possible presence of contaminants on the surfaces of the foil. Figure 2.2b 
shows a heat treated 13 ply neat Ti-15-3 laminate after creep testing at 566°C (lOOx 
magnification). The dark appearance of the specimen indicates the rather widespread 
precipitation of the alpha phase after aging. No distortion of the beta grains as a 
result of the creep test is evident. 


Ti-15-3 exhibits minimal directionality as 0.0065 inch foil [23]. The alloy has been 
shown to exhibit significant thermoviscoplastic behavior [15,37]. Final strength is 
dependent upon aging temperature and prior cold work, with lower aging 
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temperatures (to approximately 900°F minimum) and greater cold working (to 
approximately 60% maximum) producing the strongest material [36,41]. The higher 
the aging temperature, the coarser the a particles [20]. At optimum strength, 
metastable p-Ti alloys may contain up to 60% a [13]. Ti-15-3 plates less than 0.250 
inches thick can be fully solution treated by air cooling without difficulty [2,28]. 


2.3 Fibers 

The fibered laminates tested in this work had continuous SCS 6 (silicon carbide) 
fibers. The fibers were manufactured by Textron using a chemical vapor deposition 
(CVD) process. In the CVD process, a carbon monofilament fiber substrate is 
resistively heated in a tubular glass reactor. Pyrolytic graphite approximately one 
micron thick is deposited on the substrate to smooth the surface and enhance 
electrical conductivity. Then, the substrate is exposed to silane and hydrogen gases. 
The silane (silicon hydride) decomposes to form silicon carbide on the substrate, and 
the fiber is complete (see Figure 2.1) [12]. 

SCS 6 fibers are relatively oxidation resistant and retain more of their strength and 
stiffness at elevated temperatures than other types of ceramic fibers (e.g. boron). The 
room temperature strength of silicon carbide fibers ranges from 450 to 700 ksi and 
the modulus of elasticity is in the 60 Msi range [40]. Fiber properties provided by the 
manufacturer are listed in Table 2.1. 


Table 2.1 Room Temperature Physical Properties of SCS 6 fibers [12] 

Property Value 

Tensile strength 3400MPa (500ksi) 

Tensile modulus 400GPa (60msi) 


Coefficient of 
Thermal Expansion 


1.5 X10' 6 /°C (2.7x1 0" 6 /°F) 
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Figure 2.1 Photograph of "as 
unidirectional MMC . Note carbon 
fibers and silicon carbide outer 


received" Ti-15-3/SCSg 
substrate in center of 
layer. Lightest areas 


Ti-15-3. 


are 


» 



Figure 2.2a Photograph of "as received" neat Ti-15-3 
laminate etched with Kroll’s reagent (200x). The large beta 
grain structure is readily apparent. The straight lines are 
the original lamina interfaces, which were invisible prior 
to etching. 



Figure 2.2b Photograph of heat treated 13 ply neat Ti-15 
3 laminate after creep testing at 566*C (200x). The dark 
appearance indicates widespread precipitation of the alpha 
phase after aging. No distortion of the beta grains as a 
result of the creep test is evident. 






CHAPTER 3 DESCRIPTION OF THERMOVISCOPLASTIC MODELS 


3.1 General Description 

Early models of inelastic behavior were dependent upon some type of yield 
function, a plastic or inelastic flow rule, and a hardening rule, and they assumed that 
plasticity and creep could be separated. Other theories, based upon 
microphenomenological behavior or thermodynamics, have more recently been 
proposed to describe inelastic deformation. These models can be further categorized 
as either unified or uncoupled theories. The two differ in their treatment of time- 
dependent and time-independent inelastic strain components. The uncoupled 
theories separate inelastic strain into plastic and creep components; the unified 
theories do not. It should be noted, however, that there is no formal 
micromechanical basis for separating the plastic and creep components. Indeed, the 
two behaviors might be inextricably related [24]. It is this area of unified constitutive 
modeling that the present work investigates. 


3.2 Internal State Variables 


The unified constitutive models are based upon a number of state variables which 
are used to account for elastic and inelastic material behavior [6]. In the general case, 
state variables include such quantities as applied stress. Young's modulus, internal 
energy, absolute temperature, heat flux, etc., as well as several internal state variables 
(IS Vs). For the unified constitutive models in particular, the ISVs are specifically 
used to account for inelastic material behavior. Various constitutive models include 
or omit different ISV's. Three ISV's of special relevance to the models considered in 
this study are inelastic strain, back stress, and drag stress. 
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Inelastic Strain (e , ), as defined in this work, equals total strain minus elastic and 
thermal strain. It includes both plastic and creep deformation. 

Back stress (B) is a parameter used to account for kinematic hardening. During 
inelastic deformation, dislocations pile up on slip planes at barriers in the crystal. 

The pile ups produce a back stress which opposes the applied stress on the slip plane. 
When the slip direction is reversed, dislocations that previously piled up aid in 
dislocation movement. This ISV relates to kinematic hardening or the "Bauschinger 
effect" whereby the tensile plastic deformation of a metal increases the tensile yield 
strength and decreases the compressive yield strength, or vice-versa. Back stress has 
also been referred to as "rest" or "equilibrium" stress in the literature [24,27]. 


Drag stress (D) corresponds to the average dislocation density, which results in 
isotropic hardening. This ISV accounts for cyclic hardening or softening of the 
material. The effect of drag stress is similar to isotropic hardening in time- 
independent classical plasticity [15,42]. Physical mechanisms and microstructural 
features contributing to drag stress might include grain boundaries, subgrains, 
dislocation tangles, solute atoms, and precipitate particles [30]. 


The three thermoviscoplastic models considered in this study all involve two or 
more of the above ISV's; they differ only in how the evolution of the ISV's with 
time and loading is described mathematically. However, the procedures described in 
the literature for evaluating the parameters from experimental data are extremely 
vague, especially for the Miller and Walker models; a reasonably detailed (and 
usable) description was found for only the Bodner-Partom model. Therefore 
algorithms which can be used to determine the parameters involved with each 
model from experimental data were developed as a part of this study, and will be 
described in detail in this chapter. In so far as possible these algorithms have been 
verified by comparison with previous work. For example, the algorithms developed 
for this study reproduced the published Bodner-Partom parameters for Rene 95 [9] 
and In 718 [5], as well as the the published Miller and Walker parameters for 
Hastelloy at 1800°F [24], 


il 



33 The Bodner-Partom Model 


1 1 


3.3.1 Description 

The Bodner-Partom theory is a microphenomonologically based, unified 
constitutive theory for characterizing inelastic behavior. It assumes that the total 
strain rate is separable into elastic and inelastic components. The model uses two 
ISV’s: inelastic strain f|, and drag stress, D. Strain hardening is related to plastic 
work. The Bodner-Partom model assumes isotropic material behavior and 
isothermal conditions. The theory is independent of a yield criterion and loading or 
unloading conditions. There is no back stress parameter, so the theory cannot 
account for the Bauschinger effect in kinematic hardening materials. Later 
modifications to the theory include back stress as an ISV [8]. The Bodner-Partom 
model is considered especially useful for metals at elevated temperatures where the 
back stress saturates to a constant value. It is relatively simple to construct, well 
verified, and requires a minimum of experimentation. Other researchers report that 
good agreement of analytical and experimental results for commercially pure 
titanium has been obtained. 



W p = ae I 


(3.4) 



where 


e = total strain 
E| = inelastic strain 

Ej = thermal strain 

E = Young's modulus 
Wp = plastic work 
D = Drag stress 

m, Dj , D 2 , r. A, n = material constants 

Eq = limiting value of the inelastic strain rate 


The "sgn" function assigns a value of (+1) if the argument is positive and (-1) if the 
argument is negative. 


3.3.2 Data Analysis Procedures 

The Bodner-Partom model requires determination of nine constants: E, e 0 , n, m, 
A, r. Do, D], and D 2 . Two types of tests are required to determine these constants: 

• constant strain rate tensile tests at several strain rates, and 

- creep tests at several stress levels for each temperature of interest. 

Strain rates and stress levels for the above mentioned tests should be compatible 
with the values used during numerical simulation. 

A diagram of the procedure used to determine the Bodner-Partom constants is 
shown in Figure 3.1. (Figures are at the end of each chapter.) Young's modulus (E) is 
determined from the initial stress-strain response when the load is applied at test 
temperature. The constant £q is assumed to be 10^ /sec for all materials, unless the 
applied strain rates are very high. 
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The parameters n, m, D 0 , and Dj are obtained from constant strain rate tensile 
data. When stress reaches a saturated value, equation (3.2) can only be satisfied if D is 
a constant. Therefore, neglecting recovery (by rapidly loading the specimen), D can be 
assumed to be in its fully work-hardened condition and at its maximum value, Dj. 
Equation (3.2) can be written as: 



(3.5) 


Therefore, the values of n and Dj can be obtained by plotting the quanity 



versus ln(a). 


The slope of this curve equals -2n, and the y-intercept 


equals the quanity |^2n lnDj + lnj^—jj , from which the value of Dj can be 
deduced. An example of this procedure is shown in Figure 3.2. 


If work hardening is assumed to be negligible (i.e. during the initial nonlinear 
portion of the tensile response), and using rapid loading to avoid the effects of 
recovery, Eq (3.3) becomes: 

dD = m(D ] - D)dW p (3.6) 


which can be integrated to give: 


ln(D r D) = -mW p + In^ - D 0 ) 


(3.7) 


where Wp is the inelastic work and Dq is the drag stress for an annealed (i.e., non 
cold-worked) material. An expression for D as a function of stress and inelastic strain 
rate is obtained by inverting Eq (3.2): 



2n 


D = o 


(3.8) 
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Values for D can be obtained by inserting values of stress and inelastic strain rate 
(taken from several points along the nonlinear portion of the tensile response data) 
into Eq (3.8). The constants m and Do can then be obtained from Eq (3.7). A plot of 
ln(D] - D) vs Wp (from the constant strain rate tests) gives m (the slope) and Dq (the 
intercept). 

The constants A, r, and D 2 are computed from data obtained during constant stress 
creep tests. D 2 is the minimum calculated value of D, found using Eq (3.8). During 
secondary creep (i.e., where the creep rate is approximately constant), Eq (3.3) 
becomes: 



(3.9) 


(3.10) 


By plotting the quantity In [m (Dj - D) W p J versus the quantity r In j / 

obtained during constant stress creep tests, the slope and intercept define the 
constants r and A, respectively. 


3.4 The Miller Model 
3.4.1 Description 

The second model considered during this study was proposed by Miller [30, 31]. 
Like the Bodner-Partom model. Miller's constitutive model is strongly related to 
underlying microscopic physical mechanisms. Miller's theory consists of a 
hyperbolic sine strain-rate equation, which represents all of the inelastic strain, plus 
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work hardening/ recovery equations for the two history variables. It contains three 
ISV’s: inelastic strain, back stress, and drag stress. 

This model is most suitable for strongly workhardening materials, such as 
stainless steels. Miller’s model has a demonstrated ability to simulate annealing, the 
effects of warm working, cyclic hardenening, cyclic softening, the Bauschinger effect, 
and various transient phenomena [6,30]. It does not contain an explicit yield stress, 
and hence is classified as a unified theory. 

On the basis of Young and Sherby’s work [49], Miller assumes in the development 
of his model that the relationship between room temperature yield strength and 
prior creep stress is insensitive to prior creep temperature. He assumes that drag 
stress and back stress reach constant values during steady state creep. Miller also 
assumes that the steady state drag stress during creep, D ss , is independent of creep 
temperature and that D ss is related to steady state creep stress, o ss , by a simple 
relationship. These same assumptions are made in the application of Miller's model 
to the present study of Ti-15-3. 


The uniaxial differential form of Miller’s constitutive model can be summarized 
by the following four equations [24]: 


o = E(e - E] - ej) 


(3.11) 


l l n 


Ej = P c 0’ jsinh 


(Vn 


sgn (a - B) 


(3.12) 


B = H 1 E| - HTPcQ’JsinhATlBl] sgn (B) (3.13) 


D=H 2 |e,| [c 2 + |b|-^7 


- H 2 C 2 P c e'[sinh(A 2 D 3 )] n (3.14) 


where 


E = total strain 
£[ = inelastic strain 
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and 


£y = thermal strain 

E = Young's modulus 
D = Drag stress 
B = Back stress 


P c 9', m, n, At, A 2 , C 2 , Hj, H 2 , = material constants 


The parameter 0', defined by Miller in [30], is related to the activation energy for 
plastic flow. (Note: 6 ' is not calculated directly, but is instead determined as part of 
the term P c 0 '.) 

3.4.2 Data Analysis Procedure 


As discussed in [24], the material parameters for Miller's theory are found through 
a series of constant strain rate steady state hysteresis loops under fully reversed strain 
controlled conditions at a constant temperature. The procedure used to determine 
the material parameters is outlined in Figure 3.3. 


During monotonic loading under steady state conditions, the back stress and drag 
stress are both constant. Setting Eqs (3.13) and (3.14) to zero and solving for the steady 
state inelastic strain rate: 


£lss = P c B'[sinh(Ao ss )] 


n 


(3.15) 


As shown by lmbrie, Haisler, and Allen, [24] the constant A can be found from the 
relationship: 


In £4 - In £5 


'f(o 4 ) - f(o 5 )‘ 

In £j - In £5 


((Oj) - f(o 5 ) 


(3.16) 



17 


where f(x) = ln[sinh(Ax)] and Oj is the steady state asymptotic stress in a CCSR test 
conducted at strain rate £j. The value of "A” which satisfies Eq (3.16) is selected 
iteratively. Next, from Eq (3.15), the following relationship is derived: 

ln(e ss ) = In (P c 0') + n In [sinh(Ao 88 )J (3.17) 

When values of ln(e ss ) vs In [sinMAa^)] are plotted, the slope is n and the y 
intercept is In (P c 9'). 

The constant "A" does not appear explicitly in (3.12) through (3.14), but it is needed 
to calculate Aj and A 2 . As shown in [24]: 

A 2 = [A/(1-C 1 )] 3 (3.18) 


and 


A! = A/Q 


(3.19) 


where Cj is a constant obtainable from the cyclic stress vs cyclic strain data. C 2 and Hj 

are also estimated from cyclic stress vs cyclic strain data, using a trial and error best fit 
procedure. The next parameter to be determined is H 2 , which is obtained from a best 
fit of data from either constant strain rate tensile tests or creep tests. The initial 
value of back stress is assumed to be zero. The initial value of drag stress, Dq, is 
found from: 


Do 


sinh 


-l 


o ys - 0.002 Hj 

l/n-jO.677 


(E/M’) 1/n ]° 677 } 


(3.20) 
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where a„„ is the 0.2% offset yield strength at the temperature of interest and a strain 

ys 

rate e. 

Two modifications to the theory were necessary during the present study to 
improve the agreement between experiment and theory. First, a new parameter, 

was added to the initial drag stress equation (i.e., Eq 3.20): 


D n = 


H 4 (a vs - 0.002 Ht) 


sinh' 


{<e/M’> 1/n ] a677 } 


(3.21) 


Inclusion of this new parameter permitted modelling of the inverse strain rate 
sensitivity observed at 482°C. 

Second, H 4 as well as a second new parameter, H^, were added to the relationship 
for back stress rate (i.e., Eq 3.13): 

B = H 3 H 4 {Hj E[ - H]P c 0'JsinhA]|B j sgn (B)} (3.22) 

Together parameters H 3 and H 4 increase the strain rate sensitivity of the model with 
temperature, in accordance with measurements obtained during the present study at 
562 and 649°C. In order to reach the desired asymptotic stress extremes for the range 
of strain rates, the simulated CCSR curves were less "squared" then the 

experimental curves. A closer agreement between simulated and experimental CCSR 
data is obtainable by including the parameter Hg in the model. The constants and 

H , are determined from a trial and error best fit of the CCSR data. A comparison of 
4 

predicted response to experimental results is presented in Chapter 6. Thus, the 
parameters to be determined for the (modified) Miller model are E, p c 0 , n, H^, H 2 , 

Hj, H^, Aj, A 2 , C 2 , and Dq. 


K 
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3.5 The Walker Model 
3.5.1 Description 

Walker's constitutive theory 147] is based upon a nonlinear modification to a 
three-parameter solid (a spring and Voigt element in series). It contains three IS V's: 
back stress, drag stress, and inelastic strain. The growth law for the back stress 
contains both dynamic and static recovery terms. (A dynamic term differs from a 
static term in that a dynamic term is activated only in the presence of an inelastic 
strain rate.) In Walker’s theory, at "high" strain rates the static recovery term is 
assumed insignificant compared to the dynamic recovery term and the back stress 
becomes independent of strain rate. Only dynamic recovery terms are included in 
the growth law for the drag stress. Creep, relaxation and strain rate effects are 
modeled by a power law for the inelastic strain rate, and all of the material constants 
in the theory are functions of temperature and must be determined experimentally 
for each temperature of interest. Walker obtained favorable results from tests of 
nickel based Hastelloy-X at high temperatures. For uniaxial conditions the Walker 
model is represented by the following 6 equations: 


o = E(e - ej - ej) 

(3.23) 

_ |n-l 


a-B 


E l = Dn ( °- B) 

(3.24) 

. o 

B = (^ + n2)ej - (B - B - njEjJG 

(3.25) 

D = ns £] - nt> £j D - n^o (D - Dq)^ 

(3.26) 


8 T ( n 5 R \ • m * 1 

G= ^ (n 3 + n 4 R) ,n |TTry^R + 1 J + n 7 B ‘ B 


(3.27) 



dt' 


(3.28) 


R = 



8e, 


8t' 


Walker’s model requires determination of 16 constants. However, if one assumes 
drag stress to be constant (as Walker did for Hastelloy-X), the number of constants is 
reduced to nine, and the growth law for back stress (Eq 3.25) is simplified to: 


B = (nj + n2) £i - (B - B - njEjM n9 Ej + n7 


Ej) i 


|m-ll 


B-B 


(3.29) 


leaving only 9 parameters to be determined: E, n, m, nj , n 2 , n 7 , ng , B, and D. Thus, 
equations (3.23) and (3.24) remain unchanged by this simplifying assumption, while 

the drag stress is constant, and hence Eq (3.26) becomes D = 0. This simplifying 
assumption was also made during this study. 


3.5.2 Data Analysis Procedure 

The procedure used to determine the parameters for Walker’s model is outlined 
in Figure 3.4. Walker’s theory requires two types of tests: 

- fully reversed cyclic stress-strain tests over various strain rates and 
temperatures; 

- stress relaxation tests starting from different points on steady state hysteresis 
loops. 

o 

The first material parameter to be determined is B. (The reader is cautioned to 

o 

distinguish between the material parameter B and the time rate of change in back 
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stress, B.) Since the back stress saturates to a constant value as ej approaches °° 
(assuming nj = 0), Walker showed that for rapid cyclic loading: 



ft c 
^max " °max 


V ) 


(3.30) 


where a* „ and are the maximum values of tensile and compressive stress. 

max max 


For large strains. 


a = B + D(e) 1/n 


(3.31) 


In terms of steady state values of o and B in tension. 


e = 


t r* ' 
°max D max 


B 


n 


(3.32) 


If O] and 02 denote the maximum tensile stresses at the two maximum strain rates, 
E] and £ 2 , (B attains its saturated values for each strain rate at these stresses) then: 


n = 



In \ 


<*l * B, 


max 


o 2 - B, 


max 


(3.33) 


and 

D.«vO®- ,/n 


(3.34) 
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The constant B^ ax is determined using stress drop tests, as follows. A point is 

sought on the hysteresis loop such that the creep rate is initially zero when the stress 
is held at a constant value. Walker indicates that this is difficult to do because the 
creep rates are usually very small during unloading. He suggests using a relaxation 
test, in which the strain level at which the initial relaxation changes from positive 
(stress decrease) to negative (stress increase) is determined. This procedure was used 
during this study. During one CCSR test at each temperature (tests L22 through L24), 
the specimen was cyclically loaded until a steady state hysteresis loop developed. 
Then, during the unloading portion of the test, the strain was held at a 
predetermined value, allowing the specimen to relax. The stress vs time data was 
recorded to determine the initial stress rate during relaxation. After several minutes 
the specimen was again cycled until a steady state hysteresis loop developed. This 
procedure of cyclic loading followed by relaxation was repeated five times at each 
temperature, stopping at a different strain value each time. A best fit of the data from 
the five points at each temperature was subsequently used to estimate the strain at 
which the initial relaxation stress rate was zero. 

Constants n 2 and n^ are determined from the initial monotonic stress-strain 
curve. For rapid loading rates. 


n^ = 0 


and 


G — n<j R 


where R is the cumulative inelastic strain. Therefore, 


B = B + (B^ x -B) {l - exp[-n^(e - a/E)] J 


(3.35) 


Also, 


n 9 “ 



In 


B max ' B 


nL-° + D<a l/n {i 


dc j 1 /" 
Edej . 


(3.36) 
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where do/de is the instantaneous slope of the uniaxial monotonic stress-strain 

curve at the point where the stress is cr and the strain is e. Constant n 2 can then be 
estimated by: 


Bmax-B 


n 2 = 


n g 


(3.37) 


The final two constants, m and n 7 , are determined as follows. At the peaks of 

the hysteresis loops, B = o = 0 for tensile loadings, and it can be shown for these 
conditions: 


n 2 £4 = (B 4 - B) ng £4 + n 7 (B 4 - B) m (3.38) 

and 

n 2 65 = (B 5 - B) ng £5 + ny (B 5 - B) m (3.39) 


If o, and Oc denote the maximum tensile stresses at the two smallest strain rates, £4 
4 D 

and £ 5 , B 4 and B 5 (the saturated back stresses for these two strain rates) may be 
determined by equation (3.31). Then, combining (3.38) and (3.39), m is obtained by: 


Wt 

m = wl 


(3.40) 


where 


W! 


= In 


n 2 £4 - ng £4 (B 4 - B) 
— 0 
n 2 £5 - ng £5 (B 5 - B) 


(3.41) 


and 



W 2 = In 


(3.42) 


B 4 -B 

b 5 -b 


Once m is known, n 7 may be found from (3.38) or (3.39). 



Figure 3.1 
parameters . 


Determination of Bodner-Partom model 
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Figure 3.3 Determination of Miller model parameters. 
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Figure 3.4 


Determination of Walker model parameters. 







CHAPTER 4 EXPERIMENTAL PROCEDURES AND EQUIPMENT 


4.1 Test Matrix 

Tables 4.1 and 4.2 contain details of the creep, constant strain rate, and cyclic 
constant strain rate tests completed in this program. The purpose of these tests was 
to build a data base for Ti-15-3, determine the experimental parameters required for 
the three constitutive models, and to verify and evaluate each of the models for use 
in a rule of mixtures formulation for unidirectional laminates. The following 
abbreviations are used in the tables: 

SSC - single step creep 

MSC - multistep creep 

CSR - constant strain rate 

CCSR - cyclic constant strain rate 


4.2 Test Equipment and Facilities 


Testing for this program was accomplished at the University of Washington and 
at NASA Langley Research Center. In addition, a small amount of verification took 
place at the United States Air Force Academy. 


4.2.1 University of Washington 


Tests at the University of Washington were accomplished using a SATEC Model G 
creep rupture tester (#G-2860-P). The tester features a power positioning resistance 
furnace, automatic load leveling, and a choice of 6:1 or 30:1 lever arm ratio. Three 
type K Chromel-Alumel thermocouples monitored temperature in the gage length 
(center two inches) of each specimen. The temperature data was recorded on a Fluke 
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model 2200B Datalogger (#3425007). Displacements were measured by a SATEC 
Model 200 extensometer in conjunction with an ATS Linearly Varying 
Displacement Transducer (LVDT). The LVDT was used to measure displacements 
and hence strains over a 5.08 cm (2.0 in) gage length. The LVDT had a linear 
displacement range of 0.050 inches. The input voltage (6.5V) for the LVDT was 
supplied by a Fluke 3330B Programmable constant voltage calibrator (#2111011). 
Displacement data was recorded on an IBM PC/AT using an IBM PC Data 
Acquisition and Control Adaptor. 


422 NASA Langley Research Center 

Tests at NASA LARC were run using two MTS Model 810 Servo-hydraulic 
Material System Testers with Instron hydraulically activated grips. One tester 
featured an induction furnace which used two infrared pyrometers for temperature 
control, while the other had a quartz lamp heating system equipped with 
thermocouples for temperature control. Displacements were measured by quartz-rod 
extensometers having a 1 inch gage length and a ± 0.15 inch linear range. Data was 
recorded with Nicolet XF-44 Data Storage systems backed by an X-Y analog plotter. 


42 Test Specimens 
4.3.1 Creep Test Specimens 

AH of the neat specimens used for creep tests consisted of thirteen plies of Ti-15-3 
foil, hot isostatically pressed, as manufactured by Textron for NASA LARC. The 
nominal specimen dimensions for the neat specimens was 0.145 X 1.27 X 15.25 cm. 
The fibered specimens were eight plies thick with nominal specimen dimensions of 
0.18 X 1.91 X 15.25 cm. All specimen material was received in the annealed condition 
and subsequently aged at 649°C for one hour prior to testing to precipitate the alpha 
phase and stabilize the microstructure. 


i 
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43.2 Constant Strain Rate Test Specimens 

The neat specimens used for CSR tests, like those used for creep tests, consisted of 
thirteen plies of hot isostatically pressed Ti-15-3 foil. Nominal specimen dimensions 
for the neat CSR specimens were the same as described in Section 4.3.1 for creep test 
specimens. 

433 Cyclic Constant Strain Rate Test Specimens 

Because of the requirement in both Miller's and Walker's models for fully 
reversed, strain-controlled cyclic loading, the parameters for these models were not 
attainable from the thin axially loaded neat specimens which were initially available 
for this work. The thin laminated specimens, preferable because they more closely 
resemble the fibered specimens, could not be used for CCSR testing because they 
buckle when subjected to compressive loads sufficiently high to cause plastic 
deformation. Therefore, thicker specimens were required for these tests. The 
specimens used for CCSR tests were machined from a 0.64 cm thick solid (not 
laminated) plate of Ti-15-3. All CCSR specimen material was received in an 
annealed condition and subsequently aged at 649°C for one hour prior to testing. 


4.4 Test Procedures 

Before each test, specimen dimensions were recorded, test equipment was warmed 
up, and associated computer equipment for test control and data acquisition was 
programmed. Finally, in order to stabilize the microstructure and simulate use 
conditions, the test specimens were aged at 649°C for one hour. 

4.4.1 Creep Tests 

Two types of creep tests were accomplished at the University of Washington: 
single step creep (SSC) and multi-step creep (MSC). During the creep tests, constant 
tensile loads were applied to each specimen using the lever arm creep frame. The 
resulting axial strain was recorded as a function of time. SSC tests on neat specimens 



32 


were conducted at three different temperatures-482, 566, and 649°C (900, 1050, and 
1200°F)-and five different stresses~34.5, 69, 103.4, 138, and 172.4 MPa (5, 10, 15, 20 and 
25 ksi). The MSC tests were designed to represent a "complex" load history. Two 
creep tests were accomplished using fibered specimens. Details concerning the creep 
tests are given in Table 4.1. 


4.4.2 Constant Strain Rate Tests 

Two types of tests were conducted at NASA LaRC: CSR and CCSR tests. In the 
constant strain rate tests, a variable tensile load was applied such that a constant 
strain rate was induced in the specimen. The CSR tests were accomplished at five 
different strain rates (1 X lO 4 , 5 X 10" 4 , 1 X 10' 3 , 5 X 10- 3 , and 1 X 10' 2 /sec) and three 
different temperatures— 482, 566, and 649°C (900, 1050, and 1200°F)~ for a total of 15 
tests. These 15 tests were done on neat (fiberless) specimens. Additionally, CSR tests 
were completed on four fibered specimens. The neat CSR specimens were strained to 
approximately 6%, while the fibered specimens were strained to approximately 
0.67%. , 


4.4.3 Cyclic Constant Strain Rate Tests 

Fully-reversed strain-controlled loading was required for the CCSR tests. In the 
cyclic constant strain rate tests, a variable tensile load was applied such that a 
constant strain rate was induced in the specimen. At a predetermined strain limit 
(normally 0.9 or 1.0% strain, except as noted below), the load was reversed. In most 
cases, the specimens were tested for a minimum of 20 cycles, or until the hysteresis 
loop stabilized. 

Five different strain rates were used for the CCSR tests— 1 X 10"*, 5 X 10"* , 1 X 10" 3 , 5 
X 10* 3 , and 1 X 10' 2 /sec-and three different temperatures-482, 566, and 649°C (900, 
1050, and 1200°F)-for a total of 15 tests. These 15 tests were done on neat (fiberless) 
specimens. 
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As previously mentioned, three of the CCSR tests (tests L22, 23, and 24) were 
specifically modified to accomodate the Walker model. After the material response 
was stabilized during the CCSR tests, the strain was stopped at a predetermined 
value and stress relaxation behavior was measured. After a period of time the 
cycling was reinitiated until the response was again stabilized, and then the strain 
was stopped at a different value. This process was continued until relaxation was 
measured at a total of five different strain values. 

Three other CCSR tests were altered to accomodate the Miller model (tests L25, 26, 
and 27). In these tests, the strain limits were increased in five increments from ±0.2 
to ±1.2%. The specimens were cycled approximately 30 to 40 times within each of the 
limits before proceeding to the next. 


4.5 Algorithms For Application of Models 


The following relationships were used in computing the unidirectional laminate 
longitudinal (Ef) and transverse (E2) moduli: 


Ei = EfV f + E m V m 


(4.1) 


v _ E f E m 
E2 " E f V m +E m Vf 


(4.2) 


where: Ef = elastic modulus of the fibers 


E m = elastic modulus of the matrix 
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Vf = volume fraction of fibers 
V m = volume fraction of matrix 

As is typical with rule of mixtures analyses, it was assumed that the fiber and matrix 
strains were equal in the longitudinal direction, and the fiber and matrix stresses 
were equal in the transverse direction. Outlines of the computer algorithms are 
supplied below. 


4.5.1 Constant Strain Rate Tests Using Neat or 0° Unidirectional Specimens 


The algorithm used for predicting the CSR response of either neat or 0° 
unidirectional laminated specimens is shown in Figure 4.1. First, the desired theory 
for modeling the thermoviscoplastic response of the matrix is selected. As 
previously mentioned, the fibers are assumed to deform elastically. Next, the 
volume fraction of fibers (Vf) is input and the volume fraction of matrix material 
(V m ) is calculated. When simulating a test with a neat specimen, the volume 
fraction of fibers is set equal to zero. Each of the constitutive models has a number of 
constants (e.g. Young's modulus for fibers and matrix, initial value of drag stress, 
etc.), and these are input next along with program parameters. Program parameters 
for the constant strain rate tests are strain rate, time step size, and time limit. The 
longitudinal modulus (Ej) is computed using Eq 4.1, and the initial elastic strain of 
the composite (e c ), matrix stress (o m ), and fiber stress (<jf) are found. 

Once the initial steps are completed, iteration with respect to time begins. Time is 
increased incrementally as before, and for each time increment matrix incremental 
inelastic strain, £j, is calculated using the selected constitutive model and then 
summed. The matrix stress (o m ) and fiber stress (Of) are increased incrementally for 
each time step by the relations: 
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o m = o m + E m *(dt*strain rate - q) 


(4.3) 


and 


c f = Of + E f *(dt*strain rate) 


(4.4) 


The stress in the laminate is determined by a rule of mixtures combination of fiber 
and matrix stresses: 

o c - O m V„ * <J,V, <45> 

Next, the incremental changes in drag stress, back stress, and work hardening 
parameter for the selected model are computed. This process is repeated (back to {A} 
on the diagram) until the desired total strain is reached. 

4.5.2 Constant Strain Rate Tests Using 90° Unidirectional Specimens 

The algorithm used for predicting the CSR response of 90° unidirectional 
laminated specimens is shown in Figure 4.2. The initial several steps are identical to 
the algorithm described above in section 4.5.1 for the CSR response of neat or 0° 
unidirectional laminate specimens. The transverse elastic modulus (E 2 ) is calculated 

using Eq 4.2 prior to time step iteration. 


Once the initial steps are completed, time is increased incrementally as before. For 
each time increment, matrix incremental inelastic strain is calculated using the 
selected constitutive model and then summed. The incremental changes in drag 
stress, back stress, and work hardening parameter for the selected model are then 
computed. The matrix stress (a m ) and fiber stress (o f ) are assumed to be equal to the 

laminate stress which is determined by the relation: 


<j c = o c + E 2 *(dt*strain rate - V m ej) 


(4.6) 
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This process is repeated (back to 1A} on the diagram) until the desired total strain is 
reached. 


4.53 Creep Tests Using Neat or 0° Unidirectional Specimens 


A flow diagram depicting the algorithm used for predicting the creep response of 
neat or 0° unidirectional laminated specimens is shown in Figure 4.3. When 
simulating a test with a neat specimen, the volume fraction of fibers is set equal to 
zero. The initial several steps are identical to the algorithms described in sections 
4.5.1 and 4.5.2 for the CSR response of 0° and 90° unidirectional laminate specimens 
In the case of the creep test simulations, however, the algorithm requires that the 
creep stress be specified instead of the strain rate. For the case of 0° fibered laminate 
creep, the composite strain is assumed to equal the fiber and matrix strain. 


During computer simulation of creep tests using the Miller model, it was found 
that computational problems arise when the creep stress is applied instantaneously. 
In order to overcome these numerical problems in the simulations and to more 
closely model the way the creep stress is applied in real tests, a "ramp" loading 
approach is used in the simulations. That is, the constant strain rate algorithm 
(described in section 4.5.1) is used until the stress in the specimen reaches 99% of the 
desired creep stress value. At that point, the program switches over to the creep test 
algorithm (shown in Figure 4.3). 

The above steps are completed before iterating with respect to time. Time is then 
increased by an increment dt (typically from .01 to .0005 seconds, depending on the 
numerical stability of the model). For each time increment, matrix incremental 
inelastic strain, £j, is calculated using the selected constitutive model. Next, the 
matrix inelastic strain is summed and the incremental changes in drag stress, back 
stress, and work hardening parameter for the model are computed. Finally, updated 
values of o m , Of, and £ c are determined. This process is repeated (back to {A} on the 

diagram) until the desired time is reached. 


ii 


4.5.4 Creep Tests Using 90° Unidirectional Specimens 


The algorithm used for predicting the creep response of 90° unidirectional 
laminated specimens is shown in Figure 4.4. The initial several steps are identical to 
the algorithm described above in section 4.5.3 for the creep response of neat or 0° 
unidirectional laminate specimens. In the case of 90° fibered laminate creep, 
however, the initial composite elastic strain (e c ) is a rule of mixtures combination of 

fiber and matrix strains: 

e c = o*(V f /E,+ V m /E ni ) (4.7) 

The creep stress is applied using a constant strain rate approach as previously 
described in section 4.5.3. Then, time is increased incrementally, and for each time 
increment, matrix incremental inelastic strain, ej, is calculated using the selected 
constitutive model. The matrix inelastic strain is summed and a portion of the 
incremental matrix inelastic strain is added to the total laminate strain according to 
the relation: 


E total - E total + E 1 


(4.8) 


Next, the incremental changes in drag stress, back stress, and work hardening 
parameter for the selected model are computed. This process is repeated (back to {A} 
on the diagram) until the desired time is reached. 
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TABLE 4 • 1 Tests Conducted At Univ Of Washington 


Test Nominal Temp Type 


U1 

U2 

U4 

U5 

U6 

U7 

U8 

U9 

U10 

Ull 

U12 

U13 

U14 

U1 5 

U16 

U17 

U18 

U19 

U20 

U21 

U22 

U23 

U24 

U26 

U27 

U28 

U29 

U30 

1731 

1)32 

U33 

U34 

U35 

U36 

U38 

U41 


ssc 

34.5 

ssc 

34.5 

ssc 

34.5 

ssc 

34 . 5 

ssc 

34.5 

ssc 

69 

ssc 

69 

ssc 

69 

ssc 

69 

ssc 

69 

ssc 

69 

ssc 

103.4 

ssc 

103.4 

ssc 

103.4 

ssc 

103.4 

ssc 

103.4 

ssc 

103.4 

ssc 

138 

ssc 

138 

ssc 

138 

ssc 

138 

ssc 

138 

ssc 

138 

ssc 

172.4 

ssc 

172.4 

ssc 

172.4 

ssc 

172.4 

ssc 

172.4 

MSC 

48.3/96.6 

MSC 

48.3/96.6 

MSC 

48.3/96.6 

MSC 

48.3/96.6 

MSC 

48.3/96.6 

MSC 

48.3/96.6 

MSC 

♦Various 

MSC 

48.3/96.6 


482 

NEAT 

482 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

482 

NEAT 

482 

NEAT 

566 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

482 

NEAT 

482 

NEAT 

538 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

482 

NEAT 

482 

NEAT 

566 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

482 

NEAT 

566 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

482 

NEAT 

482 

NEAT 

566 

NEAT 

566 

NEAT 

649 

NEAT 

649 

NEAT 

566 

0 

566 

90 


a 


♦Test U38 had creep stresses of 48.3, 96.5, 167, and 262 
MPa. 
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Table 4.2 Tests Conducted At LARC 


TEST 

# 

lT 

L2 

L3 

L4 

L5 

L6 

L7 

L8 

L9 

L10 

LI 1 

L12 

LI 3 

L14 

L15 

L16 

LI 7 

L18 

L19 

L20 

L21 

L22 

L23 

L24 

L25 

L26 

L27 

L28 

L29 

L30 

L33 

L34 

L35 

L36 


STRAIN 

TYPE RATE ( 1 / s ) 


TEMP TYPE 


SPECIMEN 


CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 

CCSR 


1X10-4 

1X10-4 

1X10-4 

5X10-4 

5X10-4 

5X10-4 

1X10-3 

1X10-3 

1X10-3 

5X10-3 

5X10-3 

5X10-3 

1X10-2 

1X10-2 

1X10-2 

1X10-4 

1X10-4 

1X10-4 

5X10-4 

5X10-4 

5X10-4 

1X10-3 

1X10-3 

1X10-3 

5X10-3 

5X10-3 

5X10-3 

1X10-2 

1X10-2 

1X10-2 

1X10-4 

1X10-4 

1X10-4 

1X10-4 


482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

482 

566 

649 

649 

482 

566 

649 




















Choose model for matrix inelastic response 



Figure 4.3 Creep test algorithm for neat or 0' 
unidirectional specimens. 
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Figure 4.4 
specimens . 


Creep test algorithm for 90 unidirectional 










CHAPTER 5 EXPERIMENTAL RESULTS 


Much of the data collected during this study has been described in an earlier report, 
which was released in 1991 [45]. In particular, all measurements obtained using neat 
(or fiberless) specimens have been described in Ref 45. This earlier report contains 
roughly seventy plots of data which were collected during the course of this study. A 
limited number of tests were performed using fiber-reinforced specimens. The data 
collected during these tests were not presented in the earlier report and are included 
herein in Appendix A. A brief discussion of the main features associated with each 
type of test will be presented in this Chapter. 


5.1 Neat Specimen Results 

Since all of the data collected using neat specimens have been reported 
elsewhere^ 45 ', in the interests of brevity the original data plots for neat specimens 
will not be repeated in this report. Comparisons between theory and measurements 
for selected tests will be presented in Chapter 6. 

5.1.1 Creep Tests 

Data measured during the Creep tests using neat specimens are plotted in 
Appendix A, Reference 45. As described in Chapter 4 of the present report, creep tests 
were conducted at three different temperatures (482°C, 566°C, and 649°C) and five 
different stress levels ranging from 34.5 MPa to 172.4 MPa. The creep rate of the neat 
Ti-15-3 specimens increased with increasing creep stress and increasing temperature, 
as would be expected. A very rapid increase in creep rate occurred above 566°C 
and/or 103.4 MPa. The creep compliance curves imply that Ti-15-3 is linearly 
viscoplastic at temperatures at or below 482°C and nonlinear viscoplastic at the two 
higher temperatures. 

5.1.2 Constant Strain Rate Tests 

Data measured during the Constant Strain Rate Tests using neat specimens are 
plotted in Appendix B, Reference 45. During the constant strain rate tests, stress was 




45 


controlled so as to produce a predetermined strain rate. Five different strain rates 
were imposed, ranging from a low of .0001 /s to a high of .01 / sec. Once again, the 
three test temperatures of 482°C, 566°C, and 649°C were utilized in conjunction with 
each of the five strain rates. 

Interestingly, Ti-15-3 exhibited strain softening at high strain levels during all CSR 
tests; at strain levels beyond yielding, the slope of the stress-strain curve became 
negative. Also, in many of the CSR tests (as well as several of the CCSR tests) Ti-15-3 
exhibited an upper/lower yield point phenomenen. 

The initial modulus of Ti-15-3 was relatively insensitive to strain rate but 
generally decreased with increasing temperature. The modulus decreased from a 
reported value of 104.8 GPa at room temperature [26] to 70.7 GPa at 649°C, a drop of 
approximately 33 percent. The yield stress and tensile stength were found to be 
highly sensitive to both temperature and strain rate. 


5.1.3 Cyclic Constant Strain Rate Tests 

Data measured during the Cyclic Constant Strain Rate Tests using neat specimens 
are plotted in Appendix C, Reference 45. During the cyclic constant strain rate tests, 
stress was controlled so as to produce a predetermined strain rate until a limiting 
total strain was reached (usually ± 0.01 ) at which time the loading was reversed. Five 
different strain rates were imposed, ranging from a low of .0001 /s to a high of .01 /sec. 
Once again, the three test temperatures of 482°C, 566°C, and 649°C were utilized in 
conjunction with each of the five strain rates. 

As discussed in [45], the CCSR tests were extremely difficult to accomplish. Many 
of the tests had to be repeated; one was attempted five times before even marginal 
results could be attained. Additionally, during the tests in which an inductance 
furnace was used to heat the specimens the data recorded by the Nicolet XF-44 was 
severely contaminated by electronic noise. Although the source of the noise has not 
been clearly established, it is almost certainly linked to the induction furnace. 

Where possible the data obtained during these tests were smoothed using a 
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commercially available software package called Vupoint 146J. Smoothing was 
accomplished by fitting the data, segment by segment, with polynomials of order 2 
through 8. The data collected during tests in which a quartz lamp was used to 
heatthe specimens resulted in a much "quieter" electronic environment in which to 
record the CCSR data; this data did not require smoothing. 


5.2 Fibered Specimen Results 

The data collected using fiber-reinforced specimens is presented in Appendix A of 
this report. A brief description of each type of test conducted is given below. 

5.2.1 Creep Tests 

Two creep tests were completed using laminated specimens. Creep strains 
measured for a 90° laminate at 566°C and 48.3 MPa are shown in Figure A.l, and the 
creep response of a 0° laminate at 566°C and various incrementally increasing loads 
is shown in Figure A.2. Because of the high stiffness of the 0° specimen, the 
deformations experienced at low loads were too small for the extensometer and 
LVDT to accurately measure. As a result, during subsequent comparisons with 
theory a correction was made to the data by adding the expected elastic strains to the 
original data. The resulting values are plotted in Figure B.69 for a creep stress level 
of 262 MPa. 

5.2.2 Constant Strain Rate Tests 


Four constant strain rate tests were completed using laminated specimens. The 
results of L33, shown in Figure A.3, are for a 0° laminate at 649°C and a strain rate of 
.0001 m /m/sec. Figures A.4 through A.6 show the CSR response of 90° laminates at 
482, 566, and 649°C from tests L34 through L36, respectively. The stresses induced by 
the enforced strain histroy are much lower than expected-even lower than the 
corresponding stresses measured using neat specimens. The probable explanation for 
this discrepancy is that the fiber-matrix interface was extremely weak, reducing the 
stiffness of the laminates in the transverse direction. 


a 



CHAPTER 6 APPLICATION OF THE CONSTITUTIVE MODELS 


f 


6.1 Measured Model Parameters for Ti-15-3 


6.1.1 The Bodner-Partom Model 

As described in Chapter 3, the Bodner-Partom theory assumes that the inelastic 
component of deformation rate is a function of inelastic strain, £j, and drag stress, D. 
Strain hardening is related to plastic work. There is no back stress parameter in the 
form of the model used in this study, so it has been useful for metals where the back 
stress saturates to a constant value. 

The method used to determine the Bodner-Partom parameters was presented in 
detail in Chapter 3. The method requires data from constant strain rate (CSR) tests at 
several strain rates and from creep tests at several stress levels. Since the Bodner- 
Partom parameters are temperature dependent, the approach used in modelling Ti- 
15-3 was to find one set of parameters for each temperature to fit all five strain rates. 
This approach met with moderate success. It was determined that (for Ti-15-3, at 
least) the parameters are also somewhat strain-rate dependent. That is, the 
parameters that fit the data from the CSR tests having the four fastest strain rates (5 
X 10" 4 , 10" 3 , 5 X 10" 3 , and 10‘ 2 /sec) had to be adjusted to fit the data from the CSR test 
with the slowest strain rate (1 O' 4 /sec). 

Contributing to this difficulty is the fact that Ti-15-3 exhibits an upper-lower yield 
point phenomenon and strain softening. These characteristics greatly complicate the 
Bodner-Partom parameter analysis. Because the theory was under consideration for 
implementation into a "rule of mixtures" analysis for a composite laminate, the 
decision was made to use the parameters that best fit the data up to approximately 
0.8% strain. 

Table 6.1 contains the values obtained for the Bodner-Partom parameters at the 
three test temperatures and five strain rates considered. Those values marked by an 



asterisk were found to be particularily strain rate dependent; the values listed are the 
recommended "best" value for Ti-15-3 at the strain rates considered in this study. 


6.1.2 Miller Model 

Miller's constitutive model, as described in Chapters 3 and 4, requires creep, 
constant strain rate, and cyclic constant strain rate tests, including CCSR tests in 
which the strain limits are increased in increments. As expected, the Miller 
parameters for Ti-15-3 are highly temperature dependent. The material behavior is 
substantially affected by strain rate. At 900°F, Ti-15-3 shows an inverse strain rate 
sensitivity. Unlike most materials, there is more inelastic deformation at the stress 
and strain extremes at the fastest strain rate than there is at the slowest. At the two 
higher temperatures, the trend was reversed: there is more inelastic deformation at 
the stress and strain extremes at the slowest strain rates than there is at the fastest. A 
difficulty arises, however, because the asymptotic stress for each of the two highest 
temperatures at the fastest strain rate was double the stress at the slowest strain rate. 
This feature could not be accounted for using the original version of the Miller 
model, and hence the modifications described in Chapter 3 were implemented. 

Table 6.2 contains values for Ti-15-3 Miller parameters at the three test 
temperatures considered. All of the parameters were found to be rate insensitive, 
except for the parameter H 4 at 482°C. The value listedfor H 4 is the recommended 
"best" value for Ti-15-3 at the strain rates considered in this study. 


6.1.3 Walker Model 

The method used to determine the Walker parameters was described in detail in 
Chapter 3. The method requires data from cyclic constant strain rate (CCSR) tests at 
several strain rates and from stress drop or relaxation tests at several strains. As 
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expected, the Walker parameters for Ti-15-3 are highly temperature dependent. At 
566°C and 649°C the parameters were found to be essentially independent of strain 
rate, but at 482°C the parameters were somewhat strain rate dependent. 


As with the Bodner-Partom and Miller models, initial values for the Walker 
parameters were obtained using methods described above. Once the initial values 
were obtained, adjustments were made to the parameters to improve agreement 
between the experimental and predicted results. The Walker model values took less 
time to adjust than did the parameters for the other two models. Table 6.3 contains 
values for Ti-15-3 Walker parameters at the three test temperatures considered. The 
parameters were found to be strain-rate independent, except as noted by an asterisk. 


6.2 Algorithms Used To Generate CCSR, CSR, and Creep Curves 


6.2.1 Bodner-Partom Model 

The algorithm used to generate a prediction of the CCSR response of Ti-15-3 based 
on the Bodner-Partom model is shown in Figure 6.1. First, the Bodner-Partom 
model parameters are input along with certain program parameters, including strain 
rate, time step, and strain limit. Strain is then increased by an amount equal to the 
strain rate times the time increment (dt). Assuming that the initial response is 
elastic, the resulting stress is computed from Eq 3.1. 


The iterative portion of the algorithm begins at the line designated by {A}. For 
each new time increment (dt), the time rate of change of inelastic strain, plastic 
work, and drag stress are calculated. Also computed are the total inelastic strain, drag 
stress, total strain, and stress. This iterative process is continued until either the 
strain limit or the time limit is reached. Each time the strain limit is reached, the 
strain rate is reversed. When the time limit is reached, the program is ended. 
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The algorithms for generating CSR and creep curves with the Bodner-Partom 
model are identical to the CCSR algorithm with the following exceptions. For CSR 
curves, the strain rate is not reversed when the strain limit is reached. Instead, the 
program is ended. For creep tests, the CCSR approach is used until the stress is 
within one percent of the predetermined creep stress. At this point, the stress is held 
constant, and the iterative process starting at {A} in Figure 6.1 is continued. In the 
case of creep, however, the total strain increases for each time step by the 
incremental inelastic strain only; i.e., because the stress is constant, the elastic strain 
does not change as it does in the case of CSR or CCSR tests. 


6.2.2 Miller Model 

The algorithm used to generate a prediction of the CCSR response of Ti-15-3 based 
on the Miller model is shown in Figure 6.2. The process is identical to the one 
described for the Bodner-Partom model in section 6.2.1 with a few minor differences. 
First, the time rate of change of both back stress and drag stress are calculated for the 
Miller model. Second, new values for both back stress and drag stress are calculated 
at each time increment. Recall that the back stress was assumed constant in the 
Bodner-Partom model. Similarly, the algorithms for generating CSR and creep 
curves with the Miller model differ from the corresponding Bodner-Partom 
algorithms only in the requirement to calculate an incrementally changing back 
stress. 


623 Walker Model 

The algorithm used to generate a prediction of the CCSR response of Ti-15-3 is 
shown in Figure 6.3. The process is identical to the one described in section 6.2.1 for 
the Bodner-Partom model with a few minor differences. First, the drag stress is 
assumed constant with the Walker model, so the time rate of change of back stress 
and a new value for back stress are calculated at each time increment. Recall that the 
back stress was assumed constant in the Bodner-Partom model. Similarly, the 


algorithms for generating CSR and creep curves with the Walker model differ from 
the corresponding Bodner-Partom algorithms only in the requirement to calculate 
an incrementally changing back stress. 
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63 Computer Simulation of Neat Specimen Response 
63.1 Neat Specimen CSR Results 

Figures B.l through B.15 (Appendix B) compare the predictions obtained using the 
Bodner-Partom model with experimental measurement for the neat specimen CSR 
response. (The predicted CSR response using the Miller and Walker models can be 
seen by examining the CCSR results.) The Bodner-Partom model gave fairly 
accurate CSR predictions at all temperatures and strain rates with one noted 
exception: the model did not accurately predict the material's upper/lower yield 
point phenomenon exhibited under some conditions. 


63.2 Neat Specimen CCSR Results 


63.2.1 Bodner-Partom Model 


As previously mentioned, the version of the Bodner-Partom model used in this 
work does not have back stress as an ISV. Hence, the ability of the model to simulate 
cyclic loading for a material exhibiting a changing back stress is uncertain. 
Furthermore, since only CSR and creep tests are required to determine the Bodner- 
Partom constants, the data obtained from the experimental CCSR tests are not used. 
Therefore, it is of interest to assess the CCSR predictive capabilities of the Bodner- 
Partom model with regard to Ti-15-3. 


Figures B.16 through B.30 (Appendix B) compare the Bodner-Partom model's 
prediction of neat specimen CCSR response to the experimental data attained during 
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testing at NASA LaRC. The Bodner-Partom model generally did an adequate job in 
simulating the CCSR response of Ti-15-3 at the three test temperatures. The best 
agreement to experimental data seemed to be at the fastest strain rate (.01 /sec) and at 
the highest temperature (649°C). The Bodner-Partom model worked well when the 
material behaved in a nearly elastic-perfectly plastic manner. 

The largest differences between experiment and prediction came in the tests where 
the material exhibited an upper/lower yield point (e.g. Figure B.18). As mentioned 
previously, none of the models simulate this phenomenon. The ultimate goal of 
this study was to apply the three models considered to a unidirectional composite 
with ceramic fibers. Ceramic fibers generally fail at low values of strain (<1%). 
Therefore, when determining the Bodner-Partom parameters, the decision was 
made to take into account the higher value of stress associated with the upper yield 
point. As a result, the predicted CCSR stresses at the strain extremes tended to be 
greater than the experimental stresses. An alternative approach would be to fit the 
data to the lower stresses occurring at greater strains (beyond the upper yield point). 


63.2.2 Miller Model 

Figures B.31 through B.45 compare the Miller model’s prediction of neat specimen 
CCSR response to the experimental data attained during testing at NASA LaRC. 

After the modifications described in section 3.4.2 were made, the Miller model 
generally did a good job in simulating the CCSR response of Ti-15-3 at the three test 
temperatures. The best agreement to experimental data seemed to be at 482°C. At the 
two higher test temperatures, where the material behaved in a nearly elastic- 
perfectly plastic manner, the model was less accurate. This is not unexpected because, 
as mentioned in Chapter 3, the model was developed for and performs best with 
strongly workhardening materials. 
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63.2.3 Walker Model 

Figures B.46 through B.60 compare the Walker model’s prediction of neat 
specimen CCSR response to the experimental data attained during testing at NASA 
LaRC. The Walker model was less accurate in predicting the CCSR response of Ti- 
15-3 than the Miller model, but more accurate than the Bodner-Partom model. 
Predicted stresses were generally within about 10-20% of experimental values at the 
three test temperatures. Unlike the Miller model, the Walker equations were able to 
model the elastic-perfectly plastic shape characteristics of the Ti-15-3 hysteresis curve, 
but difficulty was encountered in accurately matching the stress extremes. The best 
agreement with experimental data seemed to be at the two higher test temperatures, 
where the material behaved in a nearly elastic-perfectly plastic manner. The model 
was less accurate at 482°C. 


633 Neat Specimen Creep Results 

Figures B.61 through B.66 compare the creep response predicted by the three 
models to the experimental results. Comparisons are made at two creep stresses (34.5 
and 172 MPa) and the three previously specified test temperatures. 

The models have limited accuracy in predicting the creep response of Ti-15-3. In 
each of the six cases, all of the models predict a linear creep response, representative 
of steady-state creep. Conversely, the experimental results for Ti-15-3 in all six cases 
clearly show that a non-linear creep response (sometimes referred to as "phase-one 
creep” or "primary creep”) precedes steady state creep. 


If a comparison is made on the basis of the linear region, the Walker model most 
closely matches the slope of the creep vs time curves in three of the six cases (Figures 
B.61, 63, and 66), while the Miller model matches most closely in only one (Figure 
B.65). Note that the Bodner-Partom and Walker predictions coincide on Figure B.61. 
Also, the Miller model over-predicts the steady-state creep rate in four of the six 
cases. Finally, although time dependent deformation actually occurred in all six 
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cases, the Bodner-Partom model predicts only elastic deformation in three of the six 
cases (Figures B.61, 62, and 63). 


6.4 Computer Simulation of Laminate Response 

Figures B.67 through B.71 compare experimental results to predictions obtained 
using the three models for a unidirectional laminate (Vf = .29) subjected to various 
loading conditions. In Figure B.67, the predicted 566°C creep response of a 0° 
laminate is shown for a creep stress of 262 MPa. Figure B.68 shows the 566°C 
predicted creep response of a 90° laminate at a creep stress of 48.3 MPa. Predicted 0° 
laminate CSR response at 649°C and .0001 m/m/sec is shown in Figure B.69. Figures 
B.70 through B.72 show predictions of the 90° laminate CSR response at .0001 
m/m/sec for each of the three test temperatures (482, 566, and 649°C respectively). 

All three models gave very accurate predictions of the 0° CSR results at 649°C. The 
Bodner-Partom model resulted in the most accuratepredictions for the 0° creep 
results at 262 MPa. The creep and CSR results for the 90° specimens were 
inconclusive in all cases due to the unexpected weakness of the laminate in the 
transverse direction. 
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Table 6.1 Bodner-Partom Parameters For Ti-15-3 


Parameter 

482'C 

566*C 

649‘C 

n 

9.69 

.3354 

*.212 

*>1 

m 

89.47E7 

17.27E9 

28.9E9 

2 . 357E-7 

2.012E-7 

.0001 

DO 

72.69E7 

16.28E9 

2.3E9 

r 

11.934 

*2.931 

2 .39 

a 

* . 2532 

*.031 

.2 

D 2 

E 

12 . 67E7 

3.242E9 

3.5E8 

79.2E9 

73.33E9 

72.5E9 

*a = 1.27 for 

482 ’ C , .0001/s. 

r = 6 and a - • 

55 for 

566 ' C , .0001/s 

Table 6.2 

n = .206 for 649*C, .0005 and 

Miller Parameters For Ti-15-3 

.0001/s. 

Parameter 

482 * C 

566*C 

649*C 

E (MPa) 

n 

78.0E3 

533 

75.6E6 

258 

68.5E3 

180 

°ys 

€ 

.0001 

.0001 

.0001 

n 

. 7882 

.9516 

3.73 

B c e f 

2.36E-7 

2.658E-8 

5.65E-4 

Hi 

1 . 733E5 

9.5E4 

8.685E4 

a 

.00194 

.00181 

2.7E-4 

Ho 

. 1034 

.15 

.6893 

Cl 

c 2 

.05 

.05 

.05 

.69 

.69 

.69 

«3 

1.5 

5 

6 

h 4 

*var 

1 

1 

*H a = • 6 y • 7 i 

.9, 1.1, 1.15 

for .0001, .0005, 

.001, .005, 


.01/s respectively at 482 C 
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Table 6.3 Walker Parameters For Ti-15-3 


Parameter 

482'C 

566 * C 

649*C 

E 

11 .49E6 

10.95E6 

10. 37E6 

n 

4.282 

5.577 

7.825 

m 

0 

.9878 

9.405 

no 

*29690 

1076000 

10950000 

n 7 

0 

.2146 

3.14E-30 

n 9 

*.8723 

214.3 

2028.2 

R 

2038 

-1020.5 

-258.5 

°circ 

D 0 

*770440 

154000 

69038 

* at 482 * 

C, n 2 = 296900 and 

ng = 8.723 for 

.005/s and 

.01/s; D 0 

= 539300, 462300, 

277400, and 308200 for .0005, 


.001, .005, and .01 respectively. 


it 



Figure 6.1 


Bodner-Partom Model CCSR Algorithm 
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Figure 6.3 


Walker Model CCSR Algorithm 









CHAPTER 7 CONCLUSIONS AND RECOMMENDATIONS 


7.1 General Behavior of Ti-15-3 at Elevated Temperatures 

Ti-15-3 exhibits significant inelastic deformation for creep stresses above 34 MPa at 
temperatures above 482°C. The initial modulus is dependent on temperature, 
decreasing by about 33% as temperature is increased from 20°C to 649°C. The 0.2% 
yield stress and ultimate tensile strength are both strain rate and temperature 
dependent. Strengths decrease with increasing temperature and decreasing strain 
rate. Under constant strain rate or cyclic constant strain rate conditions at elevated 
temperatures, Ti-15-3 exhibits an upper yield strength effect and strain softening at 
all temperatures tested. In the fibered laminate tests, the upper yield point 
phenomenon was observed in only the 649°C 90° CSR test (Figure B.72). 


7.2 General Assessment of the Viscoplastic Models Considered 


The mechanical response of Ti-15-3 to various loading conditions at elevated 
temperatures can be reasonably modeled by all three thermoviscoplastic constitutive 
theories considered, although obtaining the required material parameters required 
extensive efforts. In all cases the analysis was complicated by the upper yield point 
and strain softening phenomenon exhibited by Ti-15-3, which none of the models 
considered is capable of predicting. 


When combined with the simnple rule of mixtures formulation, all three models 
produced good predictions for a 0° laminate loaded under constant strain rate 
conditions as well as under constant load (creep) conditions. However, evaluation of 
the models for 90° specimens was inconclusive, since even the measured 
viscoplastic response of the 90° specimens was much higher than expected; the 
viscoplastic response measured for the 90° specimens was more pronounced than 
that measured for the neat (fiberless) specimens. It is believed that the very high 


viscoplastic response of the 90° specimens is due to a low strength fiber-matrix 
interface. The predictions obtained using all three models significantly 
underestimated the creep strain response and overestimated the constant strain rate 
stress response for the 90° specimens. These discrepancies will be further discussed 
in section 7.2.2. 1. 

7.2.1 Ease of Use 

All three models required extensive testing and considerable curve fitting. The 
parameters for each of the three models were temperature dependent, and in some 
cases strain-rate dependent. 

The Bodner-Partom model was judged easiest to use for two reasons. First, the 
procedure to determine the parameters is more clearly described in the literature 
than for the other two models. Second, the version of the Bodner-Partom model 
used does not require elevated temperature CCSR testing, which is very beneficial 
from a practical standpoint. On the other hand, the version used does not include 
back stress as an ISV, which may be a significant limitation, particularly when 
attempting to model the cyclic loading of a material in which a changing back stress 
occurs. Since the version used has no provision for changing back stress, it is most 
appropriate for materials whose back stress quickly saturates to a constant value. 
Unfortunately, this may not be apparent until testing is completed. 

Once the test matrix was completed, the parameters were determined for each 
model. Although each of the three models required a certain amount of "curve 
fitting", parameters for the Walker model were easiest to obtain for Ti-15-3. The 
nearly elastic-perfectly plastic response of Ti-1 5-3 contributed to the curve fitting 
difficulties encountered with the Miller model. 
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7.2.2 Predictive Accuracy 

7.2.2.1 Fibered Specimen Constant Strain Rate Tests 

Four unidirectional specimens were successfully tested under constant strain rate 
conditions at NASA Langley. One 0° specimen was tested at 649°C and .0001 /s (test 
L33). As shown in Figure B.69, all three constitutive models predicted the stress- 
strain response within about 10 percent, with the Miller model being slightly more 
accurate. 

Three 90° specimens were tested (tests L34, 35, and 36) under a strain rate of 
.0001 /s, one at each of the three test temperatures (482, 566, 649°C). The predicted 
response and experimental results of these tests are shown in Figures B.70 through 
B.72. In each case, the specimens were initially stiffer than neat specimens tested 
under identical conditions (as expected). However, as the applied strain increased 
the corresponding stress level did not increase as expected. Near the extremes of 
applied strain, the stresses induced were lower than that measured for neat 
specimens, and only 25-35% of the predicted values for fibered specimens. Similar 
results were obtained during creep tests using the 90° specimens (described below). 
Since the initial elastic modulus of the 90° specimens was fairly well-predicted using 
the rule of mixtures, and the 0° specimen behaved as predicted, it has been 
concluded that the pronounced viscoplastic response of the 90° specimens was due 
to a very compliant and/or weak fiber/matrix interface. This argument is 
substantiated by the results of Johnson et al [26], who report that the fiber-matrix 
interface in this MMC is very weak at room temperature. 


7.22.2 Fibered Specimen Creep Tests 

One 0° (test U38) and one 90° (test U41) specimen were creep tested at 566°C. The 
90° specimen experienced a much higher inelastic strain rate than expected, probably 
due to a weak fiber-matrix interface. As a result, inelastic strains were two to four 
times as high as predicted. Figure B.68 shows the experimental results and the 
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predicted response from each of the three models. As shown in Figure B.67, the 0° 
test and theory agreed much more closely than did the 90° test and theory. The 
Bodner-Partom model agreed most closely with experimental results, followed by 
the Miller and Walker models. All three overpredicted the inelastic deformation. 


7.3 Recommendations 


The following recommendations are submitted as a result of this study: 


1. The use of an induction furnace during the CSR and CCSR testing resulted in 
electronic noise contamination of the data and should be avoided in future testing. 
Quartz lamp heaters used in some of the CCSR tests provided a far superior 
electronic environment for high temperature testing, and is recommended. 

2. Cyclic constant strain rate testing is very difficult to perform. Two of the models 
investigated during this study (the Miller or Walker models) require CCSR test 
results to determine model paramters. It is recommended that alternate means of 
determining these parameters from experimental data, which avoid the use of CCSR 
testing, be developed. 

3. For materials in which back stress is constant, or nearly constant, then the 
Bodner Partom model in its presently used form provides accurate and easily 
obtainable results. Otherwise, the Miller and Walker models are best. From a purely 
practical standpoint, the parameters for the Walker model were easier to obtain 
during this study than those for the Miller model. Therefore, the Walker model 
would be recommended for use in modeling a material which demonstrates a 
pronounced Baushinger effect. 

4. The strength of the fiber/matrix interface in the Ti-15-3/SCS 6 MMC considered 

in this study must be improved. Once this is accomplished, additional creep and 
CSR tests, with both 0° and 90° laminates, should be performed. 
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APPENDIX A 


MEASURED VISCOPLASTIC RESPONSE OF 
UNIDIRECTIONAL SCS 6 /Ti-15-3 MMCS 



Strain (m^w) 


Creep Test U41 

48.3 HPa, 566C, 90 Degree Specimen 



Figure A. 1 Experimental creep test results of a .29 Vf 

90 laminate at 566 C and 48.3 MPa 
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Multistep Creep Test U38 
566C 



Figure A. 2 Experimental creep test results of a .29 Vf 

0 laminate at 566 C and various creep stresses. 



L33 Stress vs Strata 



Figure A.3^ 
at 649 C and 


Experimental CSR results for a 0 
0001 m/m/sec. 


laminate 






135 Stress vs Strata 



Figure A. 5 Experimental CSR results for a 90 
at 566 C and .0001 m/m/sec. 


laminate 
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L36 Stress vs Strain 
M9C 



Figure A. 6 Experimental CSR results for a 90 laminate 
at 566 C and .0001 m/m/sec. 



APPENDIX B 


COMPARISON OF PREDICTED AND 


MEASURED BEHAVIOR 


Bodner Parian Hodel Results 


4BZCi .8801 n/n/sec 



Figure B.l Bodner—Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 482 C, .0001 
m/m/ sec. 
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Bodner Ptrion Model Resalts 
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Figure B.2 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 566 C, .0001 
m/m/ sec. 
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Bodner Parian Model Results 
649C, .0001 n/u/sec 



Figure B.3 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 649 C, .0001 
m/m/sec. 
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Bodner Parian Model Results 
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Figure B.4 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 482 C, .0005 
m/m/sec . 
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Bodner Purton Model Results 
56&C, .8005 M/u/sec 



Figure B.5 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 566 C, .0005 
m/m/sec. 
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Bodner Pirtan Model Results 
48ZC, .081 u/u/sec 



Figure B.7 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 482 C, .001 
m/m/ sec. 
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Bodner Purton Model Results 
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Figure B.8 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 566 C, .001 
m/m/sec. 


566C, .001 u/u/sec 
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Bodner Parian Model Besults 


649C, ,001n/n/sec 



Figure B.9 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 649 C, .001 
m/m/ sec. 
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Bodner Firtm Model lest Its 
4B2C, .085 h/s/sbc 



Figure B.10 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 482 C, .005 
m/m/sec. 
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Bodner Partm Model lent Its 



StriiTi (m/h) 

Figure B.ll Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 566 C, .005 
m/m/sec. 
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Figure B.12 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 649 C, .005 
m/m/sec . 



Bodner Parian Model Hesults 
4B2C, .01 n/n/sec 



Strain (n/n) 

Figure B.13 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 482 C, .C 
m/m/sec. 
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Figure B.15 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CSR test at 649 C, .01 
m/m/ sec. 
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Dodner-Parttm Model CCSI Results 
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Figure B.16 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 482 C, .0001 
m/m/sec . 
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Figure B.17 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 566 C, .0001 
m/m/ sec. 
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Bodner-Fartoa Model CCSI Best Its 
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Figure B.19 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 482 C, .0005 
m/m/ sec. 
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BoJner-F&rton Model CCSI Jesuits 
566C, .B6B5n/>/sec 



Figure B.20 Bodner-Partom model prediction vs 
experimental results for Ti—15— 3, CCSR test at 566 C, .0005 
m/m/ sec. 
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BoJner-P&rioa Model CCSI lesnlts 
649Cj .8BB5n/n/sec 



Figure B.21 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 649 C, .0005 
m/m/ sec. 
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Bodner-Parton Model CCSI Results 


4B2C, .BBlm/u/sec 



Figure B.22 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 482 C, .001 
m/m/sec. 
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Bodner-P&rton Model CCSI Results 
SMC, .BBln/M/sec 



Figure B.23 Bodner-Partom model prediction vs 
experimental results for Ti-15**3, CCSR test at 566 C, .001 
m/m/sec. 
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Figure B.24 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 649 C, .001 
m/m/ sec. 
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Bodner-Ptrian CCSR Results 


566C, . 005ni/u/sec 



Figure B.26 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 566 C, .005 
m/m/sec . 
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Sodner-Ptrion CCSR Results 
M9C, .Bffijn/u/sec 



Figure B.27 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 649 C, .005 
m/m/ sec. 
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Bodner-Ptrton CCSR Results 
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Figure B.28 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 482 C, .01 
m/m/sec. 
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Bodner-Pirton CCSB Results 



Figure B.30 Bodner-Partom model prediction vs 
experimental results for Ti-15-3, CCSR test at 649 C, .01 
m/m/sec. 
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Hiller Hodel CCSR Results 
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Figure B.32 Miller model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .0001 m/m/sec. 
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Figure B.33 Miller node 
results for Ti-15-3, CCSR 
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9 C , .8801 m/u/sbc 



1 prediction vs experimental 
test at 649 C, .0001 m/m/sec 



Hiller Model CCSR Results 
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Hiller Model CCSR Results 
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Figure B.35 Miller model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .0005 m/m/sec. 
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Hiller CCSR Results 
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Figure B.40 Miller model prediction vs experimental 
results for Ti-15-3, CCSR test at 482 C, .005 m/m/sec. 
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Miller Model CCSR Results 
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Figure B.41 Miller model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .005 m/m/sec. 
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Figure B.42 Miller model 
results for Ti-15-3, CCSR t 
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Hiller CCSH Results 
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Figure B.43 Miller model prediction vs experimental 
results for Ti-15-3, CCSR test at 482 C, .01 m/m/sec. 
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Ualker CCSR Results 
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Figure B.46 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 482 C, .0001 m/m/sec. 
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Ualker CCSR Results 
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Figure B.48 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 649 C, .0001 m/m/sec. 
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Walker CCSR Results 
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Figure B.49 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 482 C, .0005 m/m/ sec. 



126 


Walker CCSR Results 
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Figure B.50 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .0005 m/m/sec. 
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Figure B.51 Walker model predictio 
results for Ti-15-3, CCSR test at 649 


l vs experimental 
C. .0005 m/m/sec 
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Walker CCSR Results 
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Figure B.56 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .005 m/m/sec. 
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Figure B.57 Walker model predictio 
results for Ti-15-3, CCSR test at 649 


i vs experimental 
C, .005 m/m/sec. 
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U&lker CCSB Results 
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Figure B.59 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 566 C, .01 m/m/sec. 
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Walker CCSR Results 
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Figure B.60 Walker model prediction vs experimental 
results for Ti-15-3, CCSR test at 649 C, .01 m/m/sec. 
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Creep Test Results and Predict ions 
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Figure B.61 Experimental results and predicted response 
for a neat specimen creep test at 482C, 34.5 MPa. 
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Creep Test Semite and Predictions 
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Figure B. 62 Experimental results and predicted response 
for a neat specimen creep test at 482C, 172 MPa. 
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Creep Test Results and Predictions 
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Figure B.63 Experimental results and predicted response 
for a neat specimen creep test at 5666C, 34.5 MPa. 
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Figure B.64 Experimental results and predicted response 
for a neat specimen creep test at 566C, 172 MPa. 
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Figure B.65 Experimental results and predicted response 
for a neat specimen creep test at 649C, 34.5 MPa. 
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Creep Test Ken Its and Predictions 
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Figure B.66 Experimental results and predicted response 
for a neat specimen creep test at 649C, 172 MPa. 
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0 Degree Lm irate Creep Test 
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Figure B.67 Experimental results and predicted creep 
response of a .29 Vf 0 laminate at 566 C and 262 MPa 
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Figure B. 68 Experimental results and pz 
response of a .29 Vf 90 laminate at 566 C 
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B Degree CSR Test 
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Figure B.69 Experimental results and predicted response 
of a .29 Vf 0 laminate at 649 C and a strain rate of .0001 
m/m/sec . 
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90 legree CSR Test 
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Figure B.70 Experimental results and predicted response 
of a .29 Vf 90 laminate at 482 C and a strain rate of 
.0001 m/m/sec. 





<MP 


147 



2.5E-3 5.0E-3 7.5E-3 

Strain (■/.) 


Figure B.71 Predicted response of a .29 Vf 90 laminate 
at 566 C and a strain rate of .0001 m/m/sec. 
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Figure B.72 Predicted response of a .29 Vf 90 laminate 
at 649 C and a strain rate of .0001 m/m/sec. 
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